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ABSTRACT 

We report the discovery of burst oscillations at the spin frequency in ten thermonuclear 
bursts from the accreting millisecond X-ray pulsar (AMXP) IGR J17511-3057. The 
burst oscillation properties are, like those from the persistent AMXPs SAX J1808.4- 
3658 and XTE J1814~338, anomalous compared to burst oscillations from intermittent 
pulsars or non-pulsing LMXBs. Like SAX J1808. 4-3658 they show frequency drifts in 
the rising phase rather than the tail. There is also evidence for harmonic content. 
Where IGR J17511-3057 is unusual compared to the other two persistent pulsars is 
that oscillations are not detected throughout all bursts. As accretion rate drops the 
bursts get brighter and their rise/decay time scales become shorter, while the oscilla- 
tion amplitude falls below the detection threshold: first in the burst peak and then also 
in the rise. None of the bursts from IGR J17511~3057 show evidence for photospheric 
radius expansion (which might be expected to suppress oscillation amplitude) which 
allow us to set an upper limit to the distance of 6.9 kpc. We discuss the implications 
of our results for models of the burst oscillation mechanism. 
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1 INTRODUCTION 

, It is thought that many low-mass X-ray binaries (LMXBs) 

■ contain rapidly rotating neutron stars. However measuring 
' the neutron star spin has been a challenge for more than 
, 30 years. So far we have been able to measure t he spin of 

■ only 25 out of more than 150 LMXBs (see, e.g.. iLiu et al.l 
' 120071 . for a catalog) thanks to episodes of X-ray pulsations 

at the neutron star spin frequency. These episodes can last 
for periods as short as a few seconds or as long as days to 
weeks. 

When pulsations occur in the persistent X-ray emis- 
sion, the sources are known as accreting millisecond X-ray 
pulsars (AMXPs) . Today there are 13 of these systems (see 
lAhamhano et all bOlOal : iMarkwardt fc StrohmaveJ |2010| . 
for the latest discovery) and their pulsations are thought 
to be powered by accretion, in a scenario in which matter 
from the accretion disk is channeled by the magnetic field 
lines onto the magnetic pole, forming a hot spot visible in 
X-rays. 
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When pulsations occur during surface thermonuclear 
explosions (known as Typ e I X-ray bursts, see, e.g 
IStrohmaver fc Bildstenl [2008, for a review), the sources 
are known as nuclear-powered X-ray pulsars (NPXPs) or 
burst oscillation sources. Prior to the discovery of IGR 
J17511-3057 there were 14 confirmed NPXPs (including the 
four which are also A MXPs, and excluding those which re- 
quire confirmation, see 
l2010bl : iGallowav era! 



Watts et ai]|2008al:lAitamhano et all 



20ld ). of which two are persistent 



AMXPs and two are intermittent AMXPs (see below). For 
clarity, in the rest of this work we will refer to the accretion- 
powered pulsations as "pulsations" and to the nuclear- 
powered pulsations as "burst oscillations". 

Burst oscillations presumably develop when the burning 
or its aftermath becomes in some way asymmetric. However 
the exact mechanism that generates burst oscillations is not 
yet known, and they are not observed in every burst even in 
sources that do show them. All theoretical models proposed 
to date fail to explain certain aspects of the data, and it 
seems increasingly likely that the magnetic field in the sur- 
face layers (largely neglected by current models) may play 
an important role. 
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In this regard systems that are both AMXPs and 
NPXPs are particularly useful. The magnetic field in these 
sources is obviously dynamically important, since it chan- 
nels accretion flow onto the stellar surface. This allows us 
to gauge field strength and any effect that channeling may 
have in establishing either te mperature or fuel distribution 
gradients in the surface layers dStrohmaver fc BildstenlbOOa 
IWatts et al.ll2008bl : IWattll2008l : iLamb et allboog fl 

Prior to the discovery of IGR J17511-3057, there 
were only four systems that were both AMXPs and 
NPXPs. If one compares the duty cycle of their pulsa- 
tions, these systems can be divided into two separate 
groups. The first group contains the AMXPs SAX J1808.4- 
3658 feiinands & van dcr Klis 19 91) and XTE J1814 - 
338 (JMarkwardt fc Swank! l2003l : IStrohmaver et all 120031 ): 
these systems s how pulsations persist e ntly during al l their 
outbursts (see iHartman et~al] l2009l : iPapitto et al] l200l 
respectively). The sec ond group c ontains the AMXPs 
HETE JI900 1-245 5 (iKaaret et al.l 2006) and Aql X-1 
IjCasella et al.l 120081 ): pulsations from these sources were 
only detected intermittently. HETE J1900. 1-2455 has been 
in outburst for almost 4.5 years, but showed pulsations 
only intermi ttently durin g the first few months of activ- 
ity (iGallowav et al.ll2007l . [2008a ). The intermittent AMXP 
Aql X-1 is particularly extreme: pulsations were detected 
only for ~ 150 sec out of the « 1.5 Msec of available data 
l|Casella et al.ll2008l ). 

These two groups of sources differ not only in their pul- 
sation duty cycle, but also in the characteristics of their 
burst oscillations. The intermittent pulsars show burst os- 
cillations with si milar chara cteristics to those observed in 
all non-AMXPs (IZhang et a l. 1998; Galloway ct al. 2008b; 
I Watts et al.ll2009al ). Burst oscillations are detected in only 
a subset of bursts. They show no harmonic content and are 
detected more often during decay than burst rise. Frequency 
drifts, which are common, are slow. Burst oscillation ampli- 
tudes tend to increase with energy. The persistent pulsars, 
by contrast, are both rather anomalous. They show burst 
oscillations in all bursts, throughout the whole burst (i.e. 
starting at the rise, continuing at the peak - except dur- 
ing photospheric-radius expansion peaks (PRE)- and dur- 
ing the decay). SAX J1808. 4-3658 shows very large fast 
frequency drifts in the burst rise t hat overshoot the spin 
frequency (|Chakrabartv et al.|[2003l ). while XTE J1814-338 
is remarkable for its fre q uency stability and phase- locking 
IjStrohmaver et al.1 l2003l : IWatts et al.1 l2005l . l2008bl l. XTE 
J1814-3 38 has burst oscillations w hose amplitude rises with 
energy (jWatts fc Strohmaverll2006r ) , and the oscillations also 
have substantial harmonic content. 

These differences support the idea that magnetic field 
plays an important part in determining burst oscilla- 
tion properties, either directly or by generating composi- 
tion/temperature/fuel depth gradients in the surface lay- 
ers. Understanding how burst oscillation properties vary be- 
tween and within the groups of persistent AMXPs, inter- 
mittent AMXPs and non-AMXPs is the first step to shed- 
ding light onto the mechanism (or mechanisms) that pro- 
duce brightness asymmetries during X-ray bursts. In this 
paper we report the discovery of burst oscillations from the 
persistent AMXP IGR J17511-3057. We present a detailed 
analysis of the burst oscillation properties and discuss how 
they fit into the existing population. 



1.1 IGR J17511-3057 

IGR J17511-3057 was first unambiguously identified by the 
INTEGRAL gamma ray mission on September 12th, 2009 
(jBaldovin et al.ll2009l ). However, the source was initially de- 
tected one day before by RXTE PGA Bulge scans as a 
rising flux attributed to two known sources in the same 
1° field of view: XTE J1751-305 and G RS 1747-312 (the 
former is a known « 435 Hz AMXP, see iMarkwardt et al.l 
[2003). The first PGA pointed observation was performed on 
Sep tember 12th, and reveal ed coherent pulsations at Ri245 
Hz (JMarkwardt et al.l2009ar ). Quasi-simultaneously with the 
PGA observation. Swift observed in the direction of the 
AMXP XTE J1751-305 and did not detect any source at 
its known Ghandra pos ition, or any other positio n in the 22 
arcmin diameter FoV (JMarkwardt et al] l2009al ) . All these 
observations together confirmed that the source in outburst 
was not the known AMXP XTE J1751-305 but a new one. 
A preliminary orbital solution yielded a pulse frequency of 
244.8337(1) Hz, an orbital period of 207.4(8) minutes, and 
a projected semi-major axis of 274(1) It-ms. The mass func- 
tion was 0.00107(2) Mq, giving a minimum companion mass 
of Ri 0.13 Mq. 

Swift observations at the beginning of the outburst de- 
tect ed thermonuclear (Type I) X-ray bursts (IBozzo et al.l 
120091) . RXTE follow-up observations led to the detection 
of burst osc i llations during these thermonuclear bursts 
(|Watts et al.l l2009bl ). Ghandra- HETG observatioii s were 
performed on September 22nd (JNowak et al.M2009l ). lead- 
ing to a refined position [a = 17^51"'0 8°.66, & = 
— 30°5 7 41 ) . Using the new coordinates, iRiggio et al.l 
([20091) and IPapitto et al.l (|2009l ) reported an improved or- 
bital ephemeris (using R XTE and XMM-New ton data, re- 
spectively). Very recently. IPapitto et al.l (|20ld ) presented a 
detailed analysis of the X-ray spectra of IGR J17511-3057 
using 70 ksec of XMM-Newton data. Near infrared follow- 
up observations iden tified the counterpa rt at a magnitude 
of K^ = 18 .( 3 ± .1 (|Torres et al.ll2009al lbl). Very recently, 
IBozzo et al.l (|2010|) reported on the Swift monitoring obser- 
vations during the outburst. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Light curves and color diagrams 

We use data from the Rossi X-ray Timing Explorer (RXTE) 
Proport ional Gounte r Array (PGA; for instrum ent informa- 
tion see lZhang et al., 1993: Jahoda et al.ll2006l ). There were 
85 pointed observations (ObsIDs: 94041-01 and 94042-01) 
containing « 2.5 to « 30 ksec of useful data per observa- 
tion. We use the 16-s time-resolution Standard 2 mode data 
to calculate X-ray colors. Hard and soft color are defined 
as the 9.7-16.0 keV / 6.0-9.7 keV and 3.5-6.0 keV / 2.0- 
3.5 keV count rate ratio, respectively, and intensity as the 
2.0-16.0 keV count rate. Type I X-ray bursts were removed, 
background was subtracted and deadtime corrections were 
made. Golors and inten sities were n ormalized by thos e of the 
Grab Nebula (see iKuulkers et al. il994; van Straaten et al.l 
|2003| . and Table 2 in I Altamirano et al.ll2008l for average col- 
ors of the Grab Nebula per PGU). 
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2.2 Energy spectra of the persistent emission 

For the PCA, we used the Standard 2 data of PCU 2, 
which was active in all observations. The background was 
estimated using the PCABACKEST V6.0 (see FTOOLS). 
We calculated the PCU 2 response matrix for each obser- 
vation using the FTOOLS routine PCARSP VlO.l. For the 
HEXTE instrument, spectra were accumulated for cluster B 
(as cluster A stopped rocking in October 2006), excluding 
the damaged detector and averaging both rocking directions 
to measure the background spectrum. Dead time corrections 
of both source and background spectra were performed us- 
ing HXTDEAD V6.0. The response matrices were created 
using HXTRSP V3.1. For both PCA and HEXTE, we fil- 
tered out data recorded during, and up to 30 minutes after 
passage through the South Atlantic Anomaly (SAA). We 
only used data when the pointing offset from the source 
was less than 0.02 degrees and the elevation of the source 
with respect to the Ea rth was greate r than 10 degrees. Us- 
ing XSPEC V11.3.2i (|Arnaudlll996l '). we fitted simultane- 
ously the PCA and HEXTE energy spectra using the 3.0- 
25.0 keV and 20.0-200.0 keV energy bands, respectively. We 
used a model consisting of a disk blackbody and a power 
law, absorbed with an equivalent Hyd r ogen column density 
^ (|Papitto et al.ll2OO9l . [2O10l ). 



2.3 Type I X-ray bursts 

We searched the Swift X-ray telescope (XRT) and the burst 
alert telescope (BAT) data for the occurrence of X-ray 
bursts. We found a total of 3 bur sts. Spectral analy sis of 
these b ursts have b e en rep orted bv lBozzo et al.l (|2010l ). We 
refer to lBozzo et al.l (|2010 l) for more information on the Swift 
data. Here we only report on the search for burst oscillations. 
We searched the Standard 1 mode PCA data (2-60 keV, 
0.125 seconds time resolution, no energy resolution) of all 
observations and found 10 Type I X-ray bursts (see Table[T]). 
In the following sections we describe the spectral and timing 
procedures separately. 



2. 3. 1 Time resolved spectroscopy of bursts detected with 
PCA 

We created energy spectra every 0.125 and 0.250 sec from 
the Event mode (E_125us_64M_0_ls) data of all the PCUs 
that were on during the burst. Given the high count rates 
during the peak of the bursts, we corrected each energy 
spectrum for dead-time using the methods suggested by the 
RXTE teairij. For each energy spectrum, we created the 
corresponding response matrix using the latest information 
available on the response of the instrument at the relevant 
times. As is common practice, we used as background the 
energy spectrum of the persistent emission taken seconds 
before each bursly. We used a black-body model to fit the 
resulting burst spectra; this choice allows us to compare our 



^ http://heasarc.gsfc.nasa.gov/docs/xte/recipes/pca_deadtinie.htmI 
^ We used 100 sec of the persistent emission to calculate the 
spectrum. However, we found no significant differences in the fits 
when the persistent-emission before or after the burst was used, 
or when using data-segments of different lengths -between 100 
and 500 seconds. 
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Figure 1. Lower panel: Intensity (2-16 keV) normalized by the 
Crab. With blue, black, red and green arrows we mark when X- 
ray bursts were detected by Swift, RXTE, Chandra and XMM- 
Newton respectively. The flare at the end of the outburst cor- 
responds to the outburs t of a nearby AMXP XTE J1751-305 
JMarkwardt et al.ll2009al ) in the PCA field of view. Upper panel: 
Percentage of the day covered by RXTE data. This estimates are 
based on Standard 2 data. 



results with previous work on burst spectra in other sources. 
With the procedure described above we are assuming that 
the X-ray spectra, after the persistent emission has been 
subtracted, are Planckian and that the observed luminosity 
of the source is: 

L = 47rc^T*i?^ 

so the unabsorbed bolometric X-ray flux is determined using 

Ftoi = aT\R/D)\ 

where a, T, R and D are the Stefan-Boltzmann constant, 
the black-body temperature, the neutron star photosphere 
radius, and the distance to the source, respectively. The ra- 
tio (R/D)^ is the normalization of the black-body model 
(bbodyrad - see XSPEC manual for details). We note that 
from the (x^) statistical point of view, the spectra of X-ray 
bursts are generally well described by black-body enrission. 
However, the emission from the neutron star and its environ- 
ment (e.g., accretion disk) is expected to be more complex 
than simple black-body emission (see, e . g.. [ van Paradiid 
1 19821 : [London et akHTosi : iKuulkers et al.ll2003i . and refer- 
ences therein) and T is a color temperature. 



2.3.2 Burst oscillations 

Burst oscillation analysis was conducted using 125 /is time 
resolution PCA event mode data, barycentered using the 
refined posit i on (a = 17''5T"08".66, S = -30°57'4l' , see 
iNowak et al.l 12009 ). We looked for data overruns which 
sometimes occur during bright bursts, but found none. We 
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Figure 2. Hardness-intensity diagram normalized by the Crab. 
Black arrows show how the spectral characteristics evolve in time. 
The blue arrow marks the time when the last X-ray burst was 
detected (in a Swift observation). Black filled circles mark the 
times when X-ray bursts were detected with RXTE. Empty circles 
mark observations where the AMXP XTE J1751-305 was also 
active. 



then extracted data in the 2-25 keV range where the bursts 
are dominant. 

Burst oscillation fractional rms amplitudes were com- 
puted using the procedure outlined in Section 2.1 of 



IWatts et al.l (J2005l 'l. We measure the 7? statistic, an un- 
binned version of the standard power spectrum. From this 
the best estimate of the pure source-only signal power in the 
absence of measurement noise, .^fignai ~ ■2's, is taken to be 
that for which the cumulative probability /„ of obtaining the 
measured power ^measured = Zm, fn[Zs : Z,„) = 0.5. Errors 
on Zs are given by the points where jn{Za '■ Zm) is 0.159 
and 0.841 (equivalent to ±1(t for a normal distribution). The 
fractional rms amplitude r is then given by 

where A'' is the total number of photons, and Nt is the num- 
ber of background photons. 

We first computed burst average fractional rms am- 
plitude, assuming a constant frequency, for the entirety of 
each burst. To reduce the contribution from the accretion- 
powered pulsations (which will contribute to the pulsed frac- 
tion if accretion continues unaffected during the burst), we 
used data from the period where the burst flux exceeds 
twice the pre-burst level. If the accretion flow and the frac- 
tional rms amplitude of the persistent pulsations remain un- 
changed during the burst, then the magnitude of the correc- 
tion is given by 

^bur^^'bur ~r '^acc.^ 'ace 



Ns 



(2) 



where A^bur and A'^acc are the number of source photons 
from the burst and accretion processes respectively, with 
Tbur and Tacc being the fractional rms amplitudes of the two 
different processes. The total number of source photons is 

Ns =iVbur + A^acc. 

The average fractional rms amplitude for each burst, 
computed in this way, is given in Table[T] We flnd signiflcant 
drifts in the frequency of the burst oscillations during many 



of the bursts: the amplitude that we report is for the con- 
stant frequency where the Z^ statistic is maximized. Fitting 
the frequency drifts would slightly increase the measured 
amplitude. Burst average fractional rms amplitudes for the 
flrst harmonic were also computed. 

To monitor the variation in frequency and amplitude 
over the course of each burst, we computed dynamical power 
spectra using overlapping 4s windows, to detect any fre- 
quency drifts. We also calculated fractional rms amplitudes 
for overlapping segments of 5000 photons (to ensure com- 
parable statistics throughout t he bursts). This is the same 
procedure as that adopted in IWatts et al.l (|2005l ). Ampli- 
tudes are only calculated for bins where the measured power 
exceeds 14 (equivalent to a 3a single trial detection). 



3 RXTE RESULTS 

3.1 The outburst evolution and the occurrence of 
X-ray bursts 

Figure [T] shows the 2-16 keV light curve of IGR J17511- 
3057. The outburst was detected on September 12th, 2009 
(MJD 55086) by INTEGRAL and by RXTE. On this day 
IGR J17511-3057 was already about 20 mCrab; its inten- 
sity increased for about 2 more days, reaching a peak lumi- 
nosity of « 27 mCrab (corresponding to a ~4.5% (d/6.9)^ 
Eddington luminosity as estimated from spectral fits in the 
2-200 keV range). After that, the intensity decreased al- 
most linearly until MJD « 55107.5, when the intensity ap- 
pears to decrease faster in time. The weakening of the source 
was interrupted by an apparent flare which started on MJD 
~ 55111 and lasted for about 3 days. This flare-like increase 
of intensity was not from IGR J 17511-3057, but due to the 
nearby AMXP XTE J1751-304 l|Markwardt et al.ll2009bh in 
the PGA fleld of view. 

In Figure [T] we mark the episodes of X-ray bursts 
as detected by R XTE (10 bursts). S wift (3 bursts, see 
IBozzo et ahlbOlOl '). Ghandra (1 b urst. seelNowak e t al. 2003) 
and XMM-Newton (2 bursts, see lPapitto et al.ll2009l . l201ol ). 
After MJD « 55105 no other bursts were detected by any 
satellite, suggesting a drop in burst rate. However, note that 
this drop in burst rate also coincides with the moment when 
the PGA sampling decreased signiflcantly (see upper panel 
of Figure [l}. In Figure [2] we plot the hardness-intensity dia- 
gram using aU PGA available data of IGR J17511-3057. At 
the beginning of the outburst decay, the hard color remains 
constant as the intensity decreases. A break occurs when the 
intensity reaches ~ 15 mGrab; the hard color then starts to 
decrease approximately linearly with intensity. 



3.2 X-ray burst spectral characteristics 

We found 10 Type I X-ray bursts in the PGA data; all of 
them were well fltted with a black-body model (with aver- 
age x^/dof between 0.8 and 1.2 during the flrst 5-10 seconds 
of the burst). All bursts showed similar temperature, radius 
and bolometric flux proflles (see Table [l]). The temperature 
and flux proflles were all single peaked, with no evidence 
for precursors. All bursts reached their maxima within 1.2 
seconds. The rise times (from 25 % to 90% of peak flux, as 
defined in iGallowav et al.ll2008br ) of all bursts were between 
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Figure 3. Temperature, radius and luminosity evolution for the 
briglitest burst in our sample (Burst 9, see Table[T]l. In the lowest 
panel we show the x^ /d.of of our fits. The Radius and Luminosity 
were calculated assuming a distance of 6.9 kpc 



0.6 to 1.2 sec. The maximum blackbody temperature {kT) 
was always between 2.5 and 3.5 keV and peak bolometric 
fluxes were in the (3 — 6.7) x 10~* erg s~^ cm~^ range. In all 
cases, the black-body radius remained approximately con- 
stant for the first ~ 3 seconds of the burst, then increased 
until reaching an average value of « 9 km (assuming a dis- 
tance of 6.9 kpc). All bursts were also similar as far as decay 
shape. Unlike other cases where double-exponential or power 
law decays are favored, a simple exponential gives a statis- 
tically good fit in most cases. Decay times were in the range 
5-8s. The total energy of each burst Etot was always in the 
range (2.5 - 3) x 10^^(d/6.9fcpc)^ erg. The burst timescales 
T, defined as the ratio between Etot / Lpeak, were all in the 
range 5.5-9s. In Figure Owe show the temperature, radius 
and luminosity profiles for the brightest burst in our sample 
(i.e., Burst 9, see Table[l]). 




9 10 11 



Persistent 2-200 Luminosity (10* erg/sec) 



Figure 4. Burst peak black body temperature (upper panel), 
bolometric luminosity at the peak (middle panel) and burst rise 
time (lower panel) as a function as of the 2-200 keV luminosity of 
the source before the occurrence of the X-ray bursts. To estimate 
the luminosity we assumed a distance of 6.9 kpc. Rise times were 
calculated using 0.1 s resolution. 



In FigureUwe show the black-body temperature (upper 
panel), the bolometric luminosity at the peak of each burst 
(middle panel) and the burst rise time (lower panel) against 
the 2-200 keV luminosity of the source before the bursts, 
respectively. Both the temperature and the bolometric lu- 
minosity appear to anti-correlate with the persistent lumi- 
nosity of the source; however, there is a clear correlation be- 
tween the burst rise time and the persistent luminosity. To 
further investigate the significance of the anti-correlations, 
we calculated the Spearman rank correlation coefficients; we 
found coefficients of —0.948 and —0.936 for the temperature 
and the burst-peak luminosity, respectively (the two-sided 
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significance of their deviations from zero equals 2.910"^ and 
6.710~^, respectively). We also studied the significance of 
the correlation between rise time and persistent luminosity 
and found Spearman rank correlation coefficient of +0.941, 
with a 5.0110"^ two-sided significance of its deviations from 
zero. Clearly, the apparent relations in Figure 3] are all sig- 
nificant. 



3.3 Distance to IGR J17511-3057 

None of the bursts showed indications of photospheric ra- 
dius expansion. By using the highest measured bolometric 
peak flux of 6.7 x 10~* erg s~^ cm~^ we can estimate an 
upper limit on the distance. We find a distance D < 6.9 kpc 
when using the empirically determined Eddington luminos- 
ity of 3.79 X 10"^* erg s~^ (this LEdd is an average value 

estimated from PRE bursts from sources situated in glob- 
■ I 1 

Kuulkers et al . 2003). This is consistent 
with the e stimations from XM M-Ncwton data recently re- 
ported bv IPapitto et ahl (|20ld ). Using the same model as 
iGallowav et al.l (2008b|jj and standard values for the mass 
and the radius of the neutron star (i.e. Mns = IAMq and 
R=\Q km), we found D < 5.76 kpc and D < 4.4 kpc for hy- 
drogen mass fractions of X = and X = 0.7, respectively. 
Higher values of the radius and the mass give higher up- 
per limits (at the very conservative case of Afjvs ~ 2.2Mq, 
71 = 15 km and X = 0, we find D < 7.1 kpc). 



3.4 Burst oscillations 

Burst oscillations are detected significantly in all bursts ob- 
served from this source, at a frequency very close to the 
known spin frequency. The properties of the burst oscilla- 
tions in IGR J17511-3057 are summarized in Table [T] In 
Figure [5] we show the fractional rms amplitude (upper pan- 
els) and the dynamical power spectra for each burst (lower 
panels) as a function of time. 

Burst oscillations are detected throughout the bursts 
only for the weaker bursts seen in the earlier part of the out- 
burst. As the intensity (and probably accretion rate) drops 
and the bursts become brighter, and the characteristics of 
the burst oscillations change. Burst oscillation amplitude 
drops below the detectability threshold first in the peak of 
the bursts, and then also in the rise. The fractional rms 
amplitude of the burst oscillations, when detected, lies be- 
tween 5 and 15%. Amplitude tends to be stronger in the 
tail than during the rise and peak. The burst-integrated 



^ The approximation was recently used by IGallowav et al.l 
(l2008a) to compare a sample of more than a thousand X-ray 
burst from different sources. The distance is given by: 



D = 8.6- ( 



Fluxp 



3-10 erg cm 



-)- 



-1/2 



( Mns u/2 
"- 1.4M« -' 



(i±^)-l/2.(l + X)-i/2kpc 



where Mpfg is the mass of the neutron star in solar masses, X is 
the mass fraction of hydrogen in the neutron star atmosphere and 
z(R) is the term that takes into account the gravitational redshift 
at the photosphere (were l + z{R) = (1 -2GMns /Rc'^)~^^^ , G is 
the gravitational constant, c the s peed of light and R the radius 
measured at the photosphere - see IGallowav et al.ll2008bl ). 



fractional amplitude (see dotted horizontal line in the top 
panels of Figure [5]) drops as the persistent luminosity de- 
creases and the X-ray bursts become stronger. This drop 
can be attributed to the fact that the oscillations are no 
longer detected throughout the burst. Harmonic content is 
detected during some bursts (see Table [l]); this is consistent 
with what has been found for the AMXI^ XTE J1814-338 
IStrohmaver et al.ll2003l : IWatts et al.]|2005l ') 

The burst oscillations show evidence for frequency drifts 
of ~ 0.1 Hz in the rise. The frequency tends to start be- 
low the spin frequency, reach a maximum in the peak and 
then stabilize very close to the spin frequency in the tail. 
In some cases the burst oscillation frequency in the tail ap- 
pears to be higher than the spin frequency (see for example 
bursts 3, 9 and 10), however this is potentially within the 
error bars on our contours and also within possible errors 
on spin frequency due to the timing noise. A detailed anal- 
ysis of the frequency drifts and the relation to the spin fre- 
quency (phase- locking), including any evidence for overshoot 
of the spin frequency (as observed for SAX J1808. 4-3658, see 
IChakrabartv et al.ll2003l ) will be presented in a companion 
paper (Patruno et al. in prep). 



4 SWIFT XRT/BAT RESULTS 

4.1 Upper limits on burst oscillations in the XRT 
data 



Of the 3 Swift bursts (see, e.g., IBozzo et al.ll2010l ). we ex- 
amined bursts 1 and 2, which occurred at UTC times 2009- 
09-14T00:50:27 and 2009-09-15T17:17:19, respectively. Swift 
XRT Windowed Timing (WT) mode are available for those 
bursts, with a time stamp resolution of 1.766 ms (Nyquist- 
limited sampling frequency of 283.126 Hz). We found no 
BAT "failed" triggers that coincide with these X-ray burst, 
i.e., only data from XRT instrument is available for these 
bursts. In the 1-7 keV band, the source peak count rate is 
~140 ct/s during both bursts, with a pre-burst background 
of «8.5 ct/s. We examined FFT power spectra with 2 s 
durations from T-10 s to T-l-50 s around the burst epochs 
listed above. No oscillations were detectable. At the pulse 
frequency, the 95% upper limit is 14% rms during the peak 
of the burst. The upper limit is less constraining in the tail 
of the burst, with typical values of 20-25%. 



4.2 Upper limits on burst oscillations in the BAT 
data 

We examined the Swift BAT data for the third burst, which 
was reported to occur at time 2009-09-30T18:31:57 (UTC, 
see JBozzo et al.ll2010l ). We searched for oscillations in the 
BAT event data, which has a time stamp resolution of 100 
fis. In a light curve with 1 s bins, we found a significantly 
enhanced BAT count rate in the time range T-3 to T+6 
compared to the reported burst time. Since BAT is typi- 
cally a background-dominated instrument, we attempted to 
maximize the signal to background ratio in order to maxi- 
mize sensitivity to oscillations. Automatic background sub- 
traction using the BAT mask weighting technique reduces 
the effective area by about 44%. Instead, we used the raw 
events and use bracketing 5 s background intervals on both 
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BURST 1 (2 PCUs) 



BURST 2 (1 PCU) 




Tme (s) 




BURST 3 (1 PCU) 



BURST 4 (1 PCU) 




Tme (s) 




BURST 5 (3 PCUs) 



BURST 6 (1 PCU) 





Tme (s) 
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Figure 5. Variation of fractional rms amplitude and frequency of burst oscillations during each burst (2-25 keV) - plot continued 
overleaf. Scales are identical in each panel for ease of comparison. Burst lightcurves are shown in each panel. Top panels: % RMS 
fractional amplitude computed for bins of 5000 photons, overlapped by 1000 photons, where the countrate is at least twice the pre-burst 
rate. Amplitudes are only shown where the power constitutes at least a Sa (single trial) significance. The dotted line shows the integrated 
amplitude computed assuming a constant frequency across the entire period (as reported in Table [T] Lower panels: Dynamical power 
spectra computed for bins of 4s, overlapped by 0.25s. Minimum power contour 20, increasing in increments of 20. The dotted line marks 
the spin frequency of IGR J17511— 3057 at the beginning of the X-ray burst. 
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Figure |5] continued . 



sides of the burst. We found that the energy range 12-18 
keV maximizes the signal to background ratio for this burst. 
Below about ~ 12 keV, the BAT energy calibration is less 
well known, and various effects like passive absorption and 
detector energy thresholds introduce further complications. 
Above 18 keV, there are still detectable burst counts, but 
the dominant background introduces more noise. Finally, we 
removed noisy detectors, disabled detectors, and detectors 
whose fractional area exposure to this source was less than 
50%. Detectors with less fractional exposure contribute pri- 
marily background. After these cuts, the background rate 
was 1806 ct/s in 13096 detectors (out of 32768 possible de- 
tectors). We constructed a single FFT power spectrum in 
the 9 s window described above and detected no significant 
oscillations. We set a 95% upper limit to oscillations at the 
pulse frequency of 10.4% rms of the net burst emission. We 
also constructed individual 2 s power spectra to search for 
shorter term variability at the pulse frequency, but did not 
detect any. 



5 DISCUSSION 

We have presented a detailed analysis of the ten X-ray bursts 
detected by RXTE from the recently discovered AMXP 
IGR J17511-3057. These bursts were recorded from the peak 
of the outburst onwards as accretion rate fell by a factor 
~ 2. The bursts become brighter and shorter as accretion 
rate falls. This behavior is consistent with the fraction of 



He in the bursts increasing as the accretion rate falls, as 
H has more time to burn to He via the hot CNO cycle 
prior to ignition (consistent with the <4% Eddington lu- 
mino sity/accretion rate we observe. See, e.g.. lBildsterJll998l . 
l200Ci ). 

Burst oscillations were detected in all of the bursts 
from IGR J17511-3057 recorded by RXTE. In this sense 
IGR J17511-3057 is like the other two persistent AMXPs 
and unlike both the intermittent AMXPs and the non- 
pulsars, where burst oscillations are detected in only a subset 
of bursts. This may be due to the fact that we observe the 
persistent AMXPs over a smaller range of accretion rates 
than the other sources; alternatively it may be associated 
with the persistent presence of a channeling magnetic field. 

The burst oscillations from IGR J17511-3057 lie within 
~ 1 Hz of the known spin frequency, confirming the al- 
ready well-established link between burst oscillation fre- 
quency and spin frequency. The oscillations also exhibit 
frequency drifts. We observe rapid drift in the burst rise, 
starting below the spin rate and ending up very close to 
it. There is some suggestion of overshoot which will be ex- 
plored more rigorously in a companion paper. The burst 
oscillation frequency stabilizes very close to the spin fre- 
quency early in the tail. This behavior is unusual compared 
to that seen in the intermittent and non-pulsars, where 
the drift is slower and often persists througho ut the tail 
(see, e.g. jMuno et al.ll2002al : IWatts et al.ll2009al ). The drifts 
seen in IGR J17511-3057 are similar to those seen in the 
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AMXP SAX J1808.4-3658 l|Chakrabartv et al.1 IJOO^ '): the 
highly stable burst oscillations seen i n XTE J1814-338 re- 
main unique ijstrohmaver et al.ll200 J : IWatts et al.ll2008bl l. 

Oscillations remain above the detectability threshold 
throughout the bursts for the weaker bursts. As the bursts 
get brighter, however, the oscillations fall below the de- 
tectability threshold - first in the peak, and eventually also in 
the rise. In many sources (including the AMXP SAX J1808- 
3658) the vanishing of oscillations in the peak is associated 
with the onset of PRE. PRE could in principle obs cure or 
shear out any asymmetry (see, e.g. ICummingj |2005| ) : alter- 
natively convection d uring PRE bu rsts may suppress the os- 
cillation mechanism l|Coopeill2008l ). However this is not the 
case for IGR J17511-3057, since none of the bursts show 
evidence for PRE. Another possibility is that there is strong 
frequency drift in the rise and peak: this would reduce the 
possibility of detection since our search assumes a constant 
frequency. To test this we searched for signals in shorter 
time windows than those used in Section 12.3.21 but found 
no evidence to support this hypothesis. Increased He frac- 
tion in the ignition column (as accretion rate falls so there 
is more time to burn H to He via the hot CNO cycle) may 
be important, but how this could affect the detectability or 
generation of oscillations in the rise and peak is not clear. 

The disappearance of oscillations in the peak and rise 
of bursts without PR E is also seen in some non-pulsars 
(JGallowav et al.ll2008bl ) but this is the first time that it has 
been seen in bursts from a persistent AMXIjj. This may 
point to a common mechanism for suppressing the devel- 
opment of oscillations, other than PRE, that operates in- 
dependently of magnetic field. However we note that PRE 
events lasting for less than 1-2 seconds would not have been 
detectable: this time interval is shorter than the periods of 
time (^ a few seconds) for which oscillations are below the 
detectability threshold. To our knowledge there is no the- 
oretical constraint on how short PRE events could be; the 
possibility that < 2 sec PRE events might suppress the burst 
oscillation mechanism for timescales of a few seconds cannot 
be excluded and should be studied in more detail. 

Burst oscillation amplitude reaches a maximum of 
about 15% rms at the fundamental frequency. This is com- 
parable to amplitudes seen in other sources (pulsars and 
non-pulsars). It is also similar to the amplitude of the 
accretion-powered pulsations from this source. In this re- 
spect IGR J17511-3057 behaves like the other persistent 
AMXPs with burst oscillations. Several of the burst oscil- 
lation trains from IGR J17511-3057 also show a detectable 
first harmonic. The only other source to show this consis- 
tently is the AMXP XTE J1814-338. Burst oscillations from 
the intermittent pulsars and the non-pulsars d o not in gen- 
eral ha ve detectable har monic content (see, e.g., Muno ct al. 
l2002bl . although see also lBhattacharvva fc Strohmaver 2005 
which discusses the existence of harmonic content in the ris- 
ing phase of bursts from the non-pulsar 4U 1636-536) . 



^ I Strohmaver et al.l (2003) found that oscillations disappeared in 
the peak of the brightest burst from XTE J1814-338 (which oc- 
curred at the lowest accretion rate). This burst showed some evi- 
dence for PRE, however the spectral evidence was not conclusive. 



6 CONCLUSIONS 

The burst oscillations of IGR J17511-3057 are, like those 
of the other two persistent AMXPs with burst oscillations, 
atypical compared to burst oscillations from the intermit- 
tent and non-pulsars. Particularly notable are the frequency 
drifts and the fact that several of the burst oscillation trains 
have detectable harmonic content. These differences between 
the AMXPs and the other sources suggest that the presence 
of a channeling magnetic field is important to the burst os- 
cillation mechanism, either directly or by establishing com- 
position or temperature gradients. 

IGR J17511-3057 has however also broadened the spec- 
trum of properties exhibited by AMXP burst oscillations. 
This AMXP is the first for which oscillation amplitudes drop 
below the detectability threshold in bright bursts, without 
the bursts exhibiting simultaneous PRE. This phenomenon 
is also seen in non-pulsars, and points to a mechanism for 
suppressing burst oscillations - other than PRE - that is less 
affected by the presence of a dynamically important mag- 
netic field. 
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Table 1. RXTE thermonuclear X-ray bursts in IGRJ17511-3057 



to 
> 

CO 

o 

o 
o 

•-0 



Bl 

-01-02 

2 



B2 

-01-06 

1 



B3 

-01-12 

1 



B4 

-01-07 

1 



B5 

-02-12 

3 



B6 

-02-200 
1 



1.72041237 2.27323586 2.61714558 6.61841723 9.60216584 
1.0 ±0.1 0.7 ±0.1 0.6 ±0.1 



ObsID" 

Num. of PCUs On 

Start time 

(MJD-l-55088.) 

Rise time (a) 1.2 ±0.1 1.1 ±0.1 0.8 ±0.1 

Decay timescale 

(s) 7.7 ±0.1 7.5 ±0.1 6.9 ±0.1 

Tau 

(a) 8.6 ±1.8 8.5 ±2.3 7.5 ± 1.9 

Peak count/rate'' ~ 5210 ~ 2600 ~ 2950 

Peak Luminosity'^ 
(10^** erg s^l (d/6.9kpc)2) 2.23 ± 0.22 2.15 ±0.36 2.61 ± 0.38 2.53 ± 0.38 2.67 ± 0.30 3.19 ± 0.45 
Total energy released^ 

(10^^ erg) 1.91 ±0.38 1.83 ± 0.47 1.95 ± 0.46 1.97 ± 0.48 2.10 ± 0.41 2.19 ± 0.50 

Duration 

(s) ~ 38 ~ 26 ~ 23 ~ 25 ~ 30 ~ 25 

Persistent L^c 
(lO^'' org s^l (d/6.9kpc)^) 1.04 ±0.03 0.97 ± 0.03 0.89 ± 0.03 0.89 ± 0.03 0.81 ± 0.03 0.86 ± 0.02 



0.32967163 
1.2 ±0.1 



7.4 ± 0.1 



7.8 ± 2.0 
~ 2880 



7.4 ± 0.1 



7.9 ± 1.7 
~ 9330 



6.7 ± 0.1 



6.9 ± 1.7 
~ 3250 



B7 

*-02-100 

1 

10.77723818 
0.6 ± 0.1 



i ± 1.7 
- 3300 



. 3.12 ± 0.46 



~ 23 
0.64 ± 0.01 



B8 
*-03-00 

4 

11.31361260 
0.8 ± 0.1 



6.6 ± 1.2 
~ 15000 



3.34 ± 0.29 



2.12 ±0.0.47 ~ 2.22 ± 0.36 



~ 30 
0.69 ± 0.02 



B9 
*-03-030 

1 

12.63290021 
0.6± 0.1 



5.5 ± 1.7 
~ 4000 

3.94 ± 0.08 

~ 2.14 ± 0.62 

~ 17 

0.57 ± 0.01 



BIO 

♦-03-04 

2 

13.28983887 
0.6± 0.1 



6.1 ± 1.4 
~ 7950 

~ 3.64 ± 0.50 

2.23 ± 0.46 

~ 23 

0.50 ± 0.01 



Osc. rms amplitude*^ 






(fundamental) 


10.8 ± 0.4 


11.2 ± 0. 


Osc. rms amplitude^ 






(1st harmonic) 


2 3+°-^ 


2 2 + °'' 


JVacc/JVb„S 


0.09 


0.09 



1.7 ± 0.6 



< 3.8 
0.08 



< 3.7 
0.08 



< 2.1 
0.08 



< 3.6 
0.06 



i.8 ± 0.6 



< 3.7 
0.07 



6.7± 0.3 

1.3 ± 0.3 
0.07 



":'*' stands for 94041-01. 

:The persistent emission has been subtracted. The count rates are not normalized by the number of PCUs on. The 2.5-24 kcV range \vas used. 
''iBolometric (blackbody) peak Luminosity assuming a distance of 6.9 kpc. 

: Frequency given is that for which the burst average amplitude is maximized assuming a constant frequency model (see text for details). 
^;Burst average fractional rms amplitude (in units of %) of the oscillations in the 2—25 keV range assuming a constant frequency (see text for details). 
^; RMS amplitudes are quoted ivhen measured power exceeded 14 (a 3<T detection). Where measured po^ver is loiver than this we quote 3fT upper limits on 
amplitude. Again these measurements assume a constant frequency model. ^: The ratio of accretion to burst photons during this period (see Equation[2l 
corrected for background.). 



+ 0.6 



5.2 



-0.5 
0.06 
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2 



3 



5.5 ± 0.4 



< 2.5 
0.06 
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1-v 
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